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Abstract. New colour distributions have been derived 
from wide field UBVRI frames for 36 northern bright el- 
liptical galaxies and a few lenticulars. The classical lin- 
ear representations of colours against logr were derived, 
with some improvements in the accuracy of the zero point 
colours and of the gradients. The radial range of signifi- 
cant measurements was enlarged both towards the galaxy 
center and towards the outskirts of each object. Thus, the 
"central colours" , integrated within a radius of 3 ", and the 
"outermost colours" averaged near the fiv — 24 surface 
brightness could also be obtained. Some typical deviations 
of colour profiles from linearity are described. Colour- 
colour relations of interest are presented. Very tight cor- 
relations are found between the U— V colour and the Mg2 
line-index, measured either at the galaxian center or at 
the effective radius. 
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1. Introduction 



elliptical and lenticulars, CD 



The "classical" data on the large scale colour distribu- 
tions of E-type galaxies relies on observations by Bender 
and Mollenhof (1987), Vigroux et al. (1988), Franx et al. 
(1989), Peletier et al. (1990), Goudfrooij et al. (1994), to 
quote only the papers discussing the 1-D profiles of colour 
against radius, as distinguished from studies of dust pat- 
terns. Most of these data were reconsidered by Michard 
(2000) (RMOO), in an attempt to collect a significant sam- 
ple of objects with a complete optical colour set, i.e. U— B, 
B— V, B— R and V— I in a coherent photometric system. 
This was adequate to confirm previous indications about 
the cause of colour gradients: these appear to be due es- 
sentially to population gradients within galaxies, with the 
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dust playing no important role, except in galaxies with 
central intense dust patterns. Such objects are rather rare 
among the Es. 

Similar to most spectral indices of stellar populations, 
the colours suffer from the well known age-metallicity de- 
generacy, and, except U— B or U— V, are not very sensitive 
to the two parameters. They are affected by dust, at least 
locally, or perhaps systematically in the central regions ac- 
cording to inferences based on a survey by Michard (1999) 
(RM99). On the other hand, they may be measured at 
lower surface brightnesses or larger radii than the line in- 
dices. They could therefore bring useful information to the 
study of fossil stellar populations, and further constraints 
upon models of the evolution of E galaxies. The present 
work aims to provide an enlarged sample of objects with 
complete colour data, extending farther in radius than in 
previous studies, and hopefully of improved accuracy. 

In Paper I, we present the usual information about 
the observations and data reduction, and part of the re- 
sults in tabular form. A larger set of results will be made 
available in electronic form. The frames, partly reduced, 
will be made available from the HYPERCAT database, 
Observatoire de Lyon. 

In Paper II, under the assumption that the observed 
colour gradients reflect abundance variations along the 
radius, metallicity gradients will be computed from the 
present data, using new colour-metallicity calibrations de- 
rived from multi-population models for E-galaxies. These 
metallicity gradients allow an estimation of central and 
mean metallicities. Statistics of galaxies included in our 
sample indicate that mean metallicities are about solar, 
in agreement with the study by Trager et al. (2000) based 
on spectral indices. 

Often used notations 

— r isophotal radius; r = (a6)^/^ for an ellipse of semi- 
axis a and h. 

— AuB colour gradient in U— B; Aub = 
d(U — B)/d(logr) and similar for other colours. 



2 



T.P. Idiart, et al.: New UBVRI colour distributions in E-type galajcies 



— dlE, boE, unE: subclassification of E galaxies as disky, 
boxy and undetermined. 

2. Observations 

The observations were performed with the 120cm Newto- 
nian telescope of the Observatoire de Haute-Provence, in 
three runs: April 1-11 2000, May 29-June 5 2000 and Jan- 
uary 18-29 2001, noted below as run 1, 2 and 3. Tables 1 
to 4 gives lists of the observed galaxies with some param- 
eters relevant on the observing conditions. A CCD target 
Tekl024 is mounted in the camera, giving a field of view 
of 11.6x11.6 'for a pixel size of 24 microns or 0.68 ". The 
relatively large field is a favorable feature of this system 
for the observation of colour distributions in nearby galax- 
ies. Less favourable are the rather poor seeing at the OHP, 
with the FWHM of star images usually in the 2-3"range, 
with values at 4 or more during periods of northern wind 
(mistral) , and also the sky illumination by ever increasing 
urban lights. 

The camera is unfortunately affected by the so-called 
"red-halo" effect. 

3. Data analysis 

3.1. Outline of the operations 

The analysis of the frames entails the following steps: 

1. The usual corrections for offset, the "flat-fielding" , and 
the interpolation of bad columns. As explained below 
we tried to improve the flat-fields by ad hoc "super- 
flats". 

2. The registration of the 5 frames in each passband to a 

common geometry, based on a set of measured coordi- 
nates for 6 to 12 stars. This is intended to simplify the 
derivation of colour maps if needed. 

3. The preparation of each frame for measurement, in- 
volving a final attempt to measure and correct residual 
large scale background trends, corrections for parasitic 
objects, a treatment against cosmic rays peaks, and 
the calibration against the available results of aperture 
photometry (see Sect. 3.2.3). 

4. Large errors in colour measurements may result from 
small differences between the widths of the PSFs of the 
two frames involved (see for instance Michard 1999, 
and the previous literature quoted therein). When the 
FWHMs of measured PSFs in the colour set for a given 
object differed by more than 10%, it was our practice 
to modify the frame PSFs and try to make them equal 
to that of the best frame of the group (see Sect. 3.2.4). 

5. The isophotal analysis of the V frame was performed 
according to Carter (1978), as implemented in the Nice 
technique described in Michard and Marchal (1994) 
(MM94). 

6. The correction procedures for the red halo effect in 
V— I, and eventually for the effects of different PSF 



Table 1. A list of the observations. The dates refer to the 
begining of the night. Two values Wq and Wi are given 
for the frame FWHM (in "): the first is the original one; 
the second was attained after the treatment tending to 
equalize the FWHM of all 5 frames in a colour set, this at 
the value of the best one. Notes: 2974 mistral, interfering 
bright star; 3115 mistral; 3193 interfering bright star; 3605 
mistral, 3607 also observed; 3608 NGC3607 interfering. 
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Table 2. A list of the observations (continued). See con- 
ventions above. Notes: 4125 clouds; 4278 NGC4283 also 
observed; 4387 mistral; 4406 (30/5) clouds; 4406 (31/5) 
mistral. 
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Table 3. A list of the observations (continued). Sec con- 
ventions above. Notes: 4551 NGC4550 also observed; 4552 
clouds. 
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Table 4. A list of the observations (continued). See con- 
ventions above. Notes: 5831 mistral and clouds. 
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U 


m654 


2400 


20.66 


2.55 


2.11 


id 


id 


B 


m655 


500 


21.28 


2.74 


2.16 


id 


id 


V 


m656 


250 


20.38 


2.07 


- 


id 


id 


R 


m657 


180 


20.08 


2.28 


2.11 


id 


id 


i 


m658 


150 


19.08 


2.33 


2.11 


5866 


02/06/00 


U 


m700 


2400 


21.16 


2.43 


1.93 


id 


id 


B 


m701 


500 


21.76 


2.26 


1.88 


id 


id 


V 


m702 


250 


20.80 


1.84 




id 


id 


R 


m703 


180 


20.36 


1.91 




id 


id 


i 


m704 


150 


19.37 


1.82 




5982 


03/06/00 


U 


m730 


2400 


20.89 


2.20 


2.00 


id 


id 


B 


m731 


500 


21.15 


2.30 


2.06 


id 


id 


V 


m732 


250 


20.46 


2.07 




id 


id 


R 


m733 


180 


20.18 


2.11 




id 


id 


i 


m734 


150 


19.13 


2.33 


2.06 



far wings in other colours, were performed (see Sect. 
3.2.5). The correction necessitates the crossed corre- 
lation of the V frame by the I PSF and conversely. 
This operation cancels out the errors in the colour dis- 
tribution induced by the red halo, but degrades the 
resolution. A correction to the calibrations performed 
before the convolutions is needed. 

7. Finally, colour measurements were performed along 
the previously found isophotal contours, and the av- 
erage isophotal colours tabulated. Our routine at this 
stage involves corrections to the adopted values of the 
sky backgrounds, in order to eliminate the obvious ef- 
fects of inaccuracies in these (see Sect. 3.2.6). 

8. The V surface brightness and colours have been col- 
lected in ad hoc files, and corrected for galactic absorp- 
tion (or reddening) and the K effect, according to the 
precepts and constants given in the Third Reference 



Table 5. Measurements of residual fluctuations in back- 
ground before and after the final "rectification". Unit: % 
of sky background. 



Colour U B V R i 

Flats-hsuperflats 1st run Oi 139 072 068 095 

id. 2nd run 1.90 0.77 0.45 0.54 0.79 

id. 3rd run 0.69 0.79 0.56 0.60 0.64 

Final treatment 1st run 0.36 0.34 0.30 0.29 0.28 

id. 2nd run 0.40 0.27 0.26 0.31 0.30 

id. 3rd run 0.28 0.28 0.24 0.20 0.32 



Catalogue of Bright Galaxies (RC3, de Vaucouleurs et 
al., 1991). The usual linear representation of colours 
against logr have been calculated in selected ranges, 
avoiding on the one hand the central regions affected 
by imperfect seeing corrections and/or by important 
dust patterns; and on the other the outer range pre- 
sumably affected by poor background corrections and 
residual noise. 

3.2. Detail of operations. 

Some important details of the above summarized proce- 
dures will now be discussed. 

3.2.1. The background 

The sky background of flat-fielded frames showed disap- 
pointingly large scale trends, specially in the U band. To 
improve upon this situation, it was tried to derive cor- 
rections by mapping the background of such frames, at 
least those which were not "filled" by a large galaxy. Such 
maps were found to be correlated, although less so than 
expected, and their average used as a "superflat". The 
quality of the background was then generally improved: if 
not, the superflats were not used. Note that, in the ob- 
serving nms of May 2000 and January 2001, a number 
of frames of "empty" fields were obtained (sometimes in 
moonlight hours) to contribute to the derivation of super- 
flats. 

A final improvement was obtained by measuring, dur- 
ing the treatment of each frame, a number of background 
patches, and subtracting a linearly interpolated map of 
these, instead of a constant. The background residual large 
scale fluctuations were often measured at the three steps 
of the procedure, that is, after the application of the flats, 
of the superflats, and after the flnal treatment. Table 5 
summarizes the results. It may be noted that the combi- 
nation of flats and superflats left large errors in our flrst 
and second runs, specially in the U colour. The flnal back- 
ground linear "rectification" allowed quite significant im- 
provements, as seen by comparing the upper and lower 
halves of the table. 
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If wc consider an E galaxy observed under the typical 
conditions of the present series (see above for a tabulation 
of sky background values), the final residual fluctuations 
quoted here represent local errors of less than 0.1 mag near 
the isophote /xs = 25. We will return later to the question 
of errors resulting from background uncertainties. 

3.2.2. Parasitic objects 

In galaxy photometry it is necessary to remove parasitic 
objects, stars and galaxies, that overlap the measured ob- 
ject. In the present work we used concurrently the follow- 
ing techniques: 

. replacement of pixels in a circle enclosing the "para- 
site" by a circle symmetric about the center of the galaxy. 

. replacement of the circle by another one chosen in a 
nearby area. 

. marking of the pixels to be discarded in such a way 
that they are later left aside in the measuring programmes. 

3.2.3. Calibrations 

The frames were calibrated by comparisons with the re- 
sults of aperture photometry. Our first choice was to use 
the UBVRI data of Poulain (1988) (PP88), and Poulain 
and Nieto (1994) (PN94) that are available for 26 objects 
of our survey, and are in Cousins's system, notably for R 
and I. In a few cases the data collected in the HYPER- 
CAT catalogues were used. These contain both primary 
data (those with an independent photometric calibration), 
and secondary data (actually calibrated with part of the 
primary data). Only primary data were used, selected ac- 
cording to our previous experience or prejudices. A few 
completely missing calibrations in R and I were replaced 
by values calculated from the tight correlations of V— R 
and V— I with B— V, derived from Poulain's data for E 
galaxies. 

Although the I filter in the camera is of Gunn's type, 
our photometry is transfered to Cousins's system through 

the calibration. It is assumed that the difference of pass- 
bands has no significant effect on colour gradients. 

3.2.4. Equalization of PSF FWHM's 

Many studies of colour distribution in galaxies are affected 

by errors resulting from the difference in the PSFs of the 
two frames involved in a colour measurement. Franx et 
al. (1989), Peletier et al. (1990), Goudfirooij et al. (1994) 
calculated the radial range where the errors due to "dif- 
ferential seeing" are larger than some accepted threshold, 
and discarded the corresponding colours. Our policy, for 
instance in RM99, was to correct for this effect by adjust- 
ing the two frames to have PSFs with a common FWHM. 

In the present study we tried to equalize the PSFs of 
the 5 frames in a given colour set. This is feasible if the 
5 frames are taken in rapid succession, so that the PSFs 



have similar widths. The frame with the narrowest PSF 
is selected, and we find by trial and error a narrow gaus- 
sian (or sum of two gaussians) which can restore another 
frame to the same quality, or rather the same FWHM, by 
deconvolution. The parameters at hand are the a of the 
gaussian and the number of iterations in the deconvolu- 
tion. In Table 1 to 4 wc give the original FWHM Wo of 
each frame, and the improved Wj after the procedure de- 
scribed here. Obviously this cannot lead to perfect results, 
and sometimes we find in our data the signature of " differ- 
ential seeing" , in the form of large colour variations, peaks 
or dips, within the seeing disk: these defects were edited 
out unless there was some good reason to suspect a gen- 
uine central colour anomaly, such as large dust patterns, 
or the jet of NGC4486. 

The reader may notice that two observations of 
NGC4406 are listed in Table 2, one of May 30, the other of 
May 31 2000. The first one was taken through fog and with 
average seeing, while for the other the "mistral" brought 
a clearer sky and very poor seeing. A special treatment 
was then applied: the central peak of the sharp images 
was "grafted" on the corresponding regions of the unsharp 
but deeper images. This explains why the Wi is so much 
narrower than the original Wo for the frames of May 31. 

3.2.5. "Red halo" and PSFs far wings. 

Before the start of this survey, the CCD camera on the 
telescope used was known to be affected by the "red 
halo" , an unfortunate property of thinned CCDs. The au- 
reoles surrounding stellar images arc obviously brighter 
and more extended in the I band than in B or V. Not only 
the red halo, but more generally the outermost wings of 
PSFs, were measured during our observing runs in 2000- 
1. The techniques and results are described in Michard 
(2001) (RMOl). The choice of appropriate star fields al- 
lowed us to extend the measurements up to a radius of 
nearly 3 ' , and down to a level of about 0.5.10~^ of the 
central peak. Due to the red halo effect, PSF wings in I 
may be a factor of 3 brighter in an extended radius range 
than the V ones. Much smaller but still significant differ- 
ences may also occur between the PSFs of various spectral 
bands, the V PSF wings always being fainter. The V PSF 
wings however, and all the others at the same time, were 
reinforced between our observing runs of spring 2000 and 
winter 2001, probably an effect of 10 months ageing of mir- 
rors coatings. The final output of the measurements are 
average "synthetic" PSFs in the format 512x512 pixels, or 
5.8x5.8 ', for each run and pass-band. 

A set of numerical experiments on model galaxies is 
also presented in RMOl, to illustrate the effects of these 
far wings on the observed surface brightness and colour 
distributions. The most striking effect occurs for the gra- 
dient Avi, which appears strongly positive, while it is 
negative according to the classical results of Bender and 
Mollenhof (1987), or Goudfrooij et al. (1994). More subtle 
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Table 6. Comparisons of average color gradients for different subsamples of E galaxies, before or after tentative 
corrections for the effects of PSF far wings according to RMOl. No corrections are applied if the average observed 
gradients are close to the adopted reference. N, number of objects. Au^b, etc. mean gradients. 



Subsample 


N 


Au-B 


a 


Ab-v 


a 


Av-R 


a 


Av-i 


a 


RMOO 


29 


-0.152 


0.048 


-0.061 


0.025 


-0.018 


0.030 


-0.053 


0.022 


2000 Observ. 


23 


-0.138 


0.037 


-0.064 


0.018 


+0.018 


0.015 


-f0.093 


0.047 


2000 Correc. 


23 


-0.116 


0.038 






-0.016 


0.013 


-0.048 


0.026 


2001 Observ. 


14 


-0.174 


0.045 


-0.080 


0.022 


-0.017 


0.013 


+0.040 


0.037 


2001 Correc. 


14 


-0.140 


0.036 










-0.062 


0.025 



17 18 19 20 21 22 23 24 



> 
OQ 



T 



• uncorrected 
o corrected 



- 1.2 



- 1.1 



„6' 
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o 
o 



- 1.0 



- 0.9 



- 0.8 



Fig. 1. Example of the "correction" of a colour profile 
through changes in the sky background constants. Abscis- 
sae: V surface brightness in mag. Ordinates: B— V. Uncor- 
rected: open circles. Corrected: filled circles. The changes 
of the background amount here to -0.25% in B and 0.10% 
in V, that is more than average (see Table 7) 



effects are found for other colours, with relatively small 
but definite changes in gradients. 

To correct for the consequences of the red halo, or other 
similar effects upon the colour distribution in the index 
CI— C2, frame CI is convolved with the PSF of frame C2 
and conversely. After this operation, the resulting images 
have been submitted to the same set of convolutions, one 
in the atmosphere plus instrument, the other in the com- 
puter: they lead therefore to correct colour distributions, 
but with a significant loss of resolution. As the convolu- 



tions attenuate the central regions of the galaxy, and much 
more so for the V frame convolved by the I PSF, the mean 
colours are biased: a correction to the calibrations per- 
formed before the convolutions is needed. This has been 
done by a comparison of simulated aperture photometry 
to the observed one, an operation also used to estimate 
the errors in calibration (see below). 

Since the extended PSFs are found with limited ac- 
curacy, it is necessary to discuss the validity of the cor- 
responding corrections obtained through crossed convolu- 
tions, the more so because of the obvious changes of the 
PSF far wings between run 1 and run 3. The mean values 
of the colour gradients for subsamples of E galaxies have 
been used for these checks, with the results of Table 6. 
For a subsample of 12 or more E galaxies the mean colour 
gradients and their dispersions cannot differ much, so that 
their values may be used as checks of the need for a cor- 
rection and its eventual success. Q The reference for these 
comparisons is the subsample in Michard (2000) (RMOO), 
mostly a rediscussion of the " classical" data by Peletier et 
al. (1990), Goudfrooij et al. (1994) and others. 

Looking at the Table 6, it is clear that the red halo 
introduces enormous errors in the V— I gradients, but that 
the corrections are remarkably successful in restoring the 
agreement of the results with the accepted reference, both 
as regards the mean values and the dispersion. The same 
may be said about the V— R gradients. The wings of the V 
PSF were strongly reinforced between our run 3 and run 
1 or 2, but much less so for the I and R PSF wings. As a 
result the red halo effect is less in V— I for the frames of 
run 3 and disappears in V— R. 

The situation is less clear for the U— B distributions. 
Our mean uncorrected gradients are in good agreement 
with the "classical" data, essentially from Peletier et al. 
(1990), as rediscussed in RMOO. On the other hand, the 



^ This assertion is questionable as pointed out by the referee. 
We therefore tested it by comparing the distributions of the 
gradients in subsamples of 15 to 25 objects, sorted out by NGC 
numbers, from the surveys of Peletier et al. 1990, Goudfrooij 
et al. 1994 and our discussion in RMOO. This procures 8 tests, 
according to the number of colours in each of the sources. It 
turns out that the subsamples are statistically coincident in 4 
cases. They differ by slightly more than the calculated errors 
in the others. 
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U PSF wings arc consistently above the V ones in all our 
runs, so that the true slopes of the U— B variations may 
be a bit smaller than the observed ones. This error might 
well be present in the classical observations. Incidentally, 
the data of Peletier et al. were obtained in U— R, and it is 
impossible to be sure that the far PSF wings of the used 
telescope were the same in both pass-bands! 

Similar remarks might be made about our B— V data. 
Since the mean measured gradient is the same for our 
data of the year 2000 and the adopted reference (and as 
a good B PSF is not available) we take as correct this set 
of results. For our data of 2001, there is evidence that the 
wings of the B PSF were slightly above those of the V one, 
so that the B— V gradients might also be biased upwards. 

It appears that a significant source of error in the mea- 
surement of the small colour gradients in E galaxies has 
hitherto been overlooked. It might he that small systematic 
errors, of the order of some 15-20%, are still present in the 
U—B or U—V gradients published here. Although such er- 
rors would not have significant astrophysical implications, 
control observations are planned. 

3.2.6. Measurements of isophotal colours 

Our procedure uses Carter's isophotal representation of 
the V frame. The successive contours at 0.1 mag intervals 
are fitted to each of the two frames to be compared, the 
mean surface brigtnesses calculated and the corresponding 
magnitudes and colours derived. A sliding mean smooth- 
ing is applied to the data for the outer contours and a 
graph of the colour against logr displayed. This might 
hopefully be linear, or nearly so, in the studied range. 
If it is not, it is our practice to introduce corrections to 
the provisional values of the sky background for one or 
both of the frames, in order to get rid, of the " breaks" in 
the colour-radius relation typically associated with a poor 
choice of the background constant. The reader is referred 
to the graphs published in Gondfrooij et al. (1994) for ex- 
amples of such features. In Fig.l we show a color profile 
with a rather important defect due to poor backgrounds, 
and its adopted correction. 

The introduction of such " aesthetic" corrections to the 
raw data might be criticized, since it assumes a regular 
behaviour of the colours at large logr. This is however a 
reasonable hypothesis: the introduced corrections remain 
small, as shown by the statistics of Table 7. It should be 
noted that the mean sky background values derived for our 
large field frames are more precise than in previous works 
based on small field frames, where the sky was not reached 
at all. The problem lies in the presence of residual large- 
scale background fiuctuations (see above): their effects are 
similar to those resulting from the poor evaluation of a 
constant background, and can be approximately corrected 
by the introduction of an ad hoc constant, or rather a set 
of constants, for the 5 frames in the colour set. 



Table 7. Corrections to provisional sky background val- 
ues applied to cancel out "breaks" in the run of colours 
aginst logr. The table gives the mean absolute values of 
the applied corrections. Unit: % of sky background. 



Colour 


U 


B 


V 


R 


I 


1st and 2nd run 
3rd run 


0.22 
0.14 


0.14 
0.07 


0.12 
0.02 


0.18 
0.08 


0.20 
0.14 


2 


10 






50 


100 




NGC 4473 



0.6 
0.5 
0.4 

0.3 

- 1.0 
0.9 
0.8 

^0.7 




10 



r (arcsec) 



50 



rto" 



-0.8 
0.7 
0.6 

HO.5 



1.4 
1.3 

-1 1.2 
1.1 



Fig. 2. Example of a set of "regular" colour profiles for 

NGC4473. In this case, the colours arc nearly linear in 
logr through the observed range. In this particular case 
a linear fit was made in the range 5-80 "to provide the 
tabulated data. 



The linear fit was finally performed on a range of log r 
selected so as to avoid the central regions affected by 
known dust patterns or possible residual seeing-induced 
errors, and the outermost regions visibly affected by devi- 
ations from the expected straight line. 

3.2.7. Errors 

The noise is not a significant source of error in this type of 
work, because averages can be performed upon thousands 

of pixels in the galaxy regions of low S/N ratio. Residual 
noise effects at large r can be easily recognized in sample 
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Fig. 3. Example of a set of colour profiles for NGC4125, a 
galaxy with a central dust pattern of importance index 3. 
In this case, the colours show a hump for r < 10 ", small 
in U— B but much larger in other colours. In this case, the 
linear fit was restricted to the range 10-80 " 
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Fig. 4. Example of a set of colour profiles for NGC3377, 
a galaxy with a central red hump in U— B or B— V but 
not in V— I. This suggests a metallicity effect. The fit was 
obtained in the range 8-80 ". 



plots of the data (see Fig. 1, 2 , 3 , 4). The main sources 
of errors lie: 

— in the calibration, giving a probable error ac, affect- 
ing equally all parts of a given object, but dependent 
upon the quality of the available aperture photometry 
in each colour. 

— in the "differential seeing'', i.e. the small differences 
in PSF between the frames in a colour set, also after 
the adjustments described above. This error ctq occurs 
only near the center of the objects, say in a diameter 
of twice the seeing FWHM (it would be a much larger 
range without the performed adjustments). There is no 
obvious reason for it to be colour dependent. 

— in the sky background residual large-scale fluctuations, 
or equivalently, errors in the sky background estimates. 
This error erg is strongly dependent upon colour and 
the studied radius, or rather the corresponding surface 
brightness. 

These three components of the total error will now be 
considered, together with other relevant topics. 

1. Errors in colours from calibration inaccuracies. 



When a number of measurements arc available for a 
given object, for instance 5 apertures in PP88 and 
PN94, a probable error of the resulting calibration is 
readily derived from the dispersion of these measured 
values about the results of simulated aperture pho- 
tometry, calculated from our data, i.e. V magnitudes, 
isophotal parameters and colours. For the preferred 
calibrations with Poulain's data, the computed error 
is often less than 0.01 in B— V or V— R but may rise 
to 0.02 in U— B or V— I. Still larger calibration errors, 
up to 0.04, have been estimated for objects with very 
scanty or uncertain aperture photometry. These prob- 
able errors apply to the zero point of the colour regres- 
sion as shown in Table 9, and to all colour data from 
the same object. 

Residual errors from PSF equalization. 
The residual errors after this step in our data treat- 
ment can be objectively ascertained by comparing the 
"central" B— R colours in our survey with the equiv- 
alent data from RM99, derived from high-resolution 
CFHT frames. For the present survey the "central 
colours" are the integrated colours within a radius of 
3 ". From RM99, Table 6, we find the colours at the 
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isophote of 1.5 ", which are likely similar. The statistics 
of the difference B-R(new)-B-R(RM99) are for 31 
objects in common: mean=0.01; a = 0.038. Assuming 
then that the errors are equal in the two surveys, the 
probable error associated with poor PSF adjustements 
is CTo = 0.027 in B— R. This source of error has no rea- 
son to vary significantly from one colour to another. 
It is independent of the error of calibration previously 
discussed. 

The comparison between the two surveys is made pos- 
sible because the same set of calibrations has been 
used, although the field of the CFHT frames was often 
not sufficient to use all calibrations apertures. A mi- 
nor part of the above differences may come from this 
source. In RM99, the B-R of NGC2768 is quoted too 
red by 0.08 and that of NGC3610 too red by 0.10. 

3. Errors from sky background inaccuracies 

Given ca and ec, the relative errors in the sky back- 
ground evaluation for frames A and C, A and C bee- 
ing a pair among UBVRI, the magnitude error in the 
colour A— C may be expressed as 5a~c = —2.5 log(l -I- 
g^^Q0.4AM^) + 2.51og(l + eclO°-^'^^^). Here A^a is 
the magnitude contrast between the object and the 
sky in colour A. Using average values for the colours 
of E-galaxies and for the sky brightnesses, A/i^i for 
any colour may be expressed in terms of A^y. Then, 
in the range of small 5a-c, the expression reduces to 
5a-c = lMh7{KAeA - ii'cec)10°-^'^''^ where Ka = 

Kv = 1 by definition; we find Ku — 3.1 and Ki = 
1.55, while Kb and Kji are slightly below 1. 
The e^,... are unknown, but it is feasible to get statis- 
tics of the linear combinations Ka^a — Kcec- Indeed, 
as explained above, we have adopted ad hoc correc- 
tions to provisional sky background values, in order to 
regularise the colour-log r relations. It is reasonable to 
assume that the errors left after these corrections are 
proportional to the adopted corrections themselves. We 
take Sa-c = V^ac = vi^A^A — Kccc)-, where 77 is a 
small constant and the Kac niay be derived from the 
statistics of the adopted corrections given in Table 7. 
In prac;tice somewhat different statistics have been cal- 
culated, to take into account the fact that our correc- 
tions for two colours are not necessarily uncorrelated. 
The constant 77, different for our observing runs of 2000 
and 2001, was chosen so as to get a system of errors 
compatible with the appearance of the data and also 
with the errors found for the slopes of the colour- log r 
relation. The finally adopted errors as from sky back- 
ground inaccuracies are given in Table 8. Note that this 
source of error is negligible for iiy = 22 or smaller. Pre- 
dicted errors are reduced for our run 3 as compared to 
the two others. 

4. Total errors 

The above estimated errors are independent and 
should be added quadratically. In the central region 



Table 8. Probable errors as associated with sky back- 
ground inaccuracies. Units: magnitude. Estimated errors 



are the same in V 


-R as 


in B- 


V. 






Colour 


Hv 


U-B 


U-V 


B-V 


V-1 


1st and 2nd run 


23 


0.023 


0.026 


0.008 


0.012 


id 


24 


0.058 


0.064 


0.020 


0.030 


id 


24.5 


0.093 


0.102 


0.032 


0.048 


3rd run 


23 


0.017 


0.019 


0.006 


0.010 


id 


24 


0.042 


0.048 


0.015 


0.026 


id 


24.5 


0.067 


0.077 


0.024 


0.041 


0.00 


0.05 




0.10 
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Fig. 5. Correlation between the colour gradients, with the 
signs changed. Abscissae: ^bv- Ordinates: ^ub- Com- 
pare with the similar diagram in RMOO, or with the origi- 
nal correlation diagram between As/j and i^uR in Peletier 
et al. (1990). The dispersion is clearly reduced here, which 
can only be attributed to an improved accuracy. 



with r < 6 ", the total error will be ctt = (c^c + ^o)^^^- 
In the mean region with 1% < jiy < 22 the total error 
equals the calibration error. Finally, in the outer range 
of nv > 23 one can use ar = {a^ + cr|)^/^. 

5. Errors from the corrections for red halo and other far 
wings effects 

To our knowledge, the crossed convolutions used to 
correct for the red halo and similar effects do not give 
rise to random errors, but rather to systematic errors 
due to inaccuracies in the adopted PSFs. Such prob- 
lems may be detected from the study of the distri- 
butions of the slopes of the colour-log r relations dis- 
cussed below. 

6. Errors in the slopes of the colour-log r relations 
These have been calculated by two complementary 
methods. On the one hand, we can look for the corre- 
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Fig. 6. Correlation between the colour gradients, with the 
signs changed.. Abscissae: Abv- Ordinates: Ay/. Here the 
comparison with similar diagrams in RMOO does not show 
much increase in the accuracy of the new data. 

lations between the slopes derived here and those from 
the literature, notably the data collected and discussed 
in RMOO. Assuming then that the errors are the same 
in both sources, we get an estimate of our slope errors, 
hopefully an upper limit. On the other hand we can 
consider the internal correlations between the slopes 
of the various colour- log r relations, specifically Aub 
and others with Abv- Figure 5 and 6 shows the cor- 
relations of the U— B and V— I colour gradients with 
that in B— V. The coefficients of correlation arc re- 
spectively 0.73 and 0.40. A weighted mean Gm4 of the 
slopes in the 4 colours has also been used as reference 
instead of Abv with analogous results. From the dis- 
persions of such correlations the slope errors can be 
estimated, if the error for the reference dB — V/dlogr 
or GmA is "guessed". The two techniques give results 
in very good agreement, the internal correlations indi- 
cating somewhat smaller errors. 

The probable errors of the slope estimates are then 
0.03 in U-B or U-V, 0.01 in B-V, 0.015 in V-R 
or B-R, 0.02 in V-I. 



4. Observational results 

The available data from the present study are: 

1. Table 9, giving for each object the linear representation 
of the colour-log r relations, i.e. the selected inner and 

outer radii of the fit, the zero point colour with its 
probable error ac and the slope. The probable errors 
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Fig. 7. Colour-colour diagram of U— B against B— V: the 
colours are calculated from the linear representations of 
Table 9 at the effective radius r^/S (dots), at re/2 (circles) 
and at the outermost range near /xy = 24 (stars) from 
Table 11. Although the various symbols refer to widely 
different regions of the galaxies they define a common re- 
lation. 
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Table 11. 
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for the slopes are given above. This table allows easy 
calculation of the colours at any radius, notably the 
effective radius and others as considered below. 

2. In electronic form only, the tables giving for each 
galaxy, as a function of the radius r, the V magnitude 

and the colours U— B, B— V, V— R, V— I. These ta- 
bles are presently available with the U— B indices as 
observed, and consistent with Table 9. They will be 
later made available after correction for suspected ef- 
fects of PSF far wings (see discussion in 3.2.5). 

3. Series of graphs from the above tables, showing the 
colours as a function of logr or of ^v- Examples of 
these graphs are shown here to illustrate a number of 
properties of the radius-colour relations. 

4. Table 10 gives the "central" colours according to sev- 
eral definitions: we have considered the colours inte- 
grated inside the area of radius r = 3 " , and the colours 
calculated for r = 1.5 "from the linear representations 
of Table 9. They should be nearly equal if these repre- 
sentations remain valid at small r, which is not always 
the case: see below for a description of typical devia- 
tions. 

5. Table 11 collects colours measured at the outermost 
range of the available data expressed in V magni- 
tude. This gross limit varies between /iy = 23 (for 
NGC4472) and 24.5. It is controlled by the size of the 
object and the "cleanliness" of the nearby field. 

The SAO galaxies NGC3115, 3607, 4550 and 5866 have 
been observed with the E-type sample. The corresponding 
results are given in the Tables, but they have been discarded 
from the discussion. 

4-.1. Statistical comparison with previous work 

— Assuming that colours result only from population 
variations, "perfect" correlations between zero point 
colours in different pass-bands are expected. We have 
compared colour-colour correlations for previous sur- 
veys and the present one. For B— I against B— V we find 
from Goudfrooij et al. (1994) a coefficient of correla- 
tion p = 0.61 (41 objects) compared to p = 0.94 from 
Table 9 (37 objects). For U— R against B— R the data 
of Peletier et al. 1990 (38 objects) give p = 0.58, while 
the present data leads to p = 0.94. The improvement 
is less if we compare our results with the discussion in 
RMOO where the calibration of the available photome- 
try was reconsidered. 

— The correlation between colour gradients should also 
be as good as allowed by errors of measurements. 
Again we compare the correlation coefficients between 
gradients from the literature and our results for the 
same colours. From Goudfrooij et al. (1994), consider- 
ing the gradient Abi against Abv, we find p — 0.78 
which is quite good. The result derived from Table 9 
is still better, i.e. p = 0.85. The improvement is more 



pronounced when the Aun against Agn correlation 
is calculated from Peletier et al. (1990): one obtains 
p — 0.36 only, instead of p = 0.88 with our data. Q 
— Our colour measurements could be obtained at much 
lower surface brightness (or larger radii) than in pre- 
vious work. The present data extend to 23.2 < /iy < 
25.2, with a median value near pv = 24.5 in all colours. 
The tables by Goudfrooij et al. (1994), available from 
CDS Strasbourg, are mostly hmited to pv — 22.75 
in B— V and pv = 22.25 in V— I (median values), 
due to the small fields of the frames, notably in the 
I band (and probably an ad hoc cut off). Their pub- 
lished graphs generally extend to the same radius in 
both colours. The tables from Peletier et al. (1990), 
again at the CDS, extend to pv = 23.2 in B— R and 
Pv = 22.6 in U— R (mean values) due to a systenatic 
cut-off at 10% of the sky. The printed tables may be 
still more severely truncated. 

In summary, our colour data generally extends 1.5 to 2 
magnitudes deeper than in previous works, so that " ex- 
ternal colours" refering to the level /iy = 24.5 when- 
ever possible, are presented in Table 11 with realistic 
error estimates. 

4-2. Description of isophotal colour profiles. 

Most of the profiles relating colour to logr, or equivalently 
to the surface brightness pvi are regular, meaning that 
they deviate very little from a straight line in the range of 
abscissae relevant to the present data, roughly r > 2 "and 
Pv < 24.5. In this case, the "central colours" integrated 
within r < 3 ", and the colours calculated from the linear 
representation at the average radius of r = 1.5, differ very 
little. Central colour according to these two definitions are 
given in Table 10: compare columns (5) and (6) for U— V, 
and (3) with (7) for V— R. Figure 2 gives an example of 
a set of regular colour profiles for NGC4473. Of course 
the linear colour-log r relation is only approximate and 
breaks down at small r for all galaxies observed at high 
resolution. CaroUo et al. (1997) obtained V—i maps of a 
number of E-galaxies from HST frames, disclosing minute 
colour structures in some cases. In RM99, the CFHT reso- 
lution proved sufficient for the detection and classification 
of B— R central " red peaks". These are smoothed out, at 
least partly, with the OHP seeing. 

Non regular profiles have been observed in the follow- 
ing cases: 

1. When an important dust pattern occurs near the cen- 
ter of an object, it produces a central red hump in the 
colour profile. This is the case for galaxies with the 
value of 3 for the dust pattern importance index (DPII) 
introduced in RM99, such as NGC2768, 4125 and 4374, 

The coefficients of correlation given here are of course dif- 
ferent from the previously calculated ones in Sect. 3.2.7.6 , 
which refered to other colours with smaller gradients. 
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but also for 5813 and 5831. The dust ring of NGC3607, 
type SAO, produces noticeable bumps in its colour pro- 
files. Fig. 3 gives an example of the colour profiles of 
such a centrally dusty galaxy, i.e. NGC4125. In such 
cases, the integrated colours within r = 3 "are redder 
than the extrapolated colours at r = 1.5: this "extra 
reddening" is smaller in U— B than in other colours. 
The consideration of the central reddening in V— R, a 
colour nearly insensitive to age-metallicity variations 
but sensitive to dust, gives a possibility to correct other 
colours for dust effects. This has been done, in some 
applications, for galaxies of DPII3. 

2. A few objects show a central red hump in some of 
the colour profiles, specially U— B and also eventually 
B-V. This is the case of NGC3377, 3379, 3610(?) and 
4494. Fig.4 illustrates the case of NGC3377. For such 
objects the "extra reddening" defined above is near 
zero in V— R and V— I. It is therefore permitted to 
attribute it to a metallicity increase rather than to 
dust. 

3. On the contrary, the central colours U— B, B— V may 

be bluer than the extrapolation of the linear portion of 
the profile. This occurs for objects with larger than av- 
erage colour gradients. The best such case is NGC4636; 
others are 4283, 4478 and the SAO 4550. An "extra 
blueing" in U— V is then apparent from the compar- 
ison of indices in columns (5)-(6) of Table 10. This 
is also marginally the case for such giant galaxies as 
NGC4406, 4472, 4649: according to RM99 the colour 
profile of such objects tend to flatten out near the cen- 
ter. 

4. NGC4486 is remarkable in showing a central "blue 
deep", the U— B colour beeing 0.51 at r = as com- 
pared to 0.72 near r = 6.5 ". This central feature is 
probably somehow related to the famous non-thermal 
jet of this galaxy. The jet is of course very conspicuous 
in U— B , with a peak colour near -0.1. Needless to say, 
both the central "blue deep" and the jet are affected 
by seeing (and our attempts to correct its effects). 



4-3. Correlations of interest. 

4.3.1. Correlations between colour gradients. 

The correlations between the various colour gradients have 

already been noted as useful tools to evaluate the probable 
errors in gradients. These correlations are displayed in Fig. 
5 and 6. 

The cocfhcicnt of correlation between Aub and Abv 
is 0.75; that between Ayi and Abv is only 0.40. For Ayr 
it falls to 0.20, because the errors are of the same order 
of magnitude as the V— R gradients. Imposing regression 
lines running through the origin, the relative slopes are 
< Aub/ Abv = 2.2 >, < Avi/Abv = 0.8 > and < 
Ayr/ Abv = 0.25 >. These relative slopes arc in good 
agreement with the results in RMOO, and its conclusion. 
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Fig. 10. Correlation between the U— V colour and the 
Mg2 index at the effective radius. This index, and the 
Mgb one merged in the data, are taken from Kobayashi 
and Arimoto, 1999. 



i.e. the negligible influence of dust upon colour gradient, 
is confirmed. 

The distribution of colour gradients for E-galaxies 

may be of interest. The following parameters are found: 
< Aub >= -0.152 with a = 0.044; < Abv >= -0.070 
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with <T = 0.021; < Ay^ >= -0.018 with cr = 0.012; 
< Avi >= -0.054 with a = 0.023. The dispersions are 
not much larger than the errors estimated above. The 
distributions arc asymmetric: there arc 4 objects with 
Abv larger than +1.8(7 above the mean, but none at less 
than the same deviation. These galaxies, with Abv clearly 
steeper than average, by about twice the estimated prob- 
able error of measurement, are NGC4283, 4478, 4564 and 
4636, seemingly a random collection. 

Remark: An attempt to sort the E-galaxics by flatten- 
ing as measured in MM94, and to look for some relation to 
the gradients, lead to negative results. Similarly no signifi- 
cant difference was found between diE and other galaxies. 

The SAO NGC4550 has quite exceptional gradients in 
all colours, and an admixture of dust and relatively young 
stars could be invoked to explain its properties. This ob- 
ject also has very remarkable kinematics, as first described 
by Rubin et al. (1992); a model has been proposed by Rix 
et al. (1992). The few other SOs in the present sample are 
similar to Es with regard to their colour gradients. 

4.3.2. Colour-colour diagrams 

Many colour-colour diagrams can be built from the present 

data. The indices used may be calciilated from the linear 
representations in Table 9 at the effective radius re, at a 
near center position rg/B, and at an intermediate position 
re/2. The system of here used is an average of esti- 
mates in MM94, Prugniel and Heraudeau, 1998, and the 
RC3. It is satisfactory that the various indices (B — V)^ , 
(B — V)j, (B ~ V)r /8' define a common diagram with 
the corresponding U— B. This may suggest that a common 
physical variable controls the variations inside an object, 
and the object to object changes, of the two colours. Such 
graphs readily show larger than average calibrations er- 
rors: for instance, the U— B colour of NGC3610 is clearly 
too red for its B— V. 

We have also traced colour-colour diagrams for the 
"central" colours, (Table 10), i.e. integrated in the radius 
r < 3 ". They are similar to those traced with the in- 
terpolated colours, but with larger dispersions: this is not 
surprising since the central colours suffer from larger er- 
rors (see 3.2.7). 

Finally, one can trace colour-colour diagrams with the 
"outermost colours" collected in Table 11. They are in fair 
agreement with the diagrams derived from Table 9, and 
extend these towards the blue. We show in Fig. 7 a com- 
posite colour-colour diagram U— B against B— V, using the 
colours ar re/8, re/2 and the outermost range. Similarly 
Fig. 8 displays the diagram of V— I against U— V. 

4.3.3. The U— V colour as a metallicity index for 
ellipticals ? 

Burstein et al. (1988) showed a correlation between the 
central Mg2 index and a global B— V colour measured in 



a large aperture (see also Bender et al. , 1993). This type 
of correlation is reconsidered here using the U— V colour 
which is much more sensitive to metallicity than B— V, and 
taking advantage of recent estimates of the Mg2 index far 
from center. 

Two correlations between the Mg2 index and U— V are 
considered in Fig. 9 and 10. The first shows the relation 
between the two quantities near the galaxy center: the 
Mg2c index is taken from the tabulation by Faber et al. 
1989. The (U — V)^ index is the integrated colour in a 
circle of radius r = 3 ". The value for three galaxies with 
important central dust patterns have been corrected by 
reference to the central bump in (V— R), a colour sensitive 
to dust but less so to metallicity changes. The coefficient of 
correlation reaches 0.825. Taking (U — V),^ as x and Mg2c 
as y we find the regression y = 0.221 ±. 027.x -0.061 ±.003. 

The second, i.e. Fig. 10, displays the correlation be- 
tween the U— V colour and the Mg2 index at the ef- 
fective radius r^. Mg2re has been taken from Kobayashi 
and Arimoto (1999) (KA99). To increase the number of 
data points the Mgi, gradients were introduced, using 
the linear relation Mgi, = \bMg2 derived from the dis- 
tribution of the values of both indices in the tables of 
KA99. The coefficient of correlation still reaches 0.72. 
Again, with the colour in x and the Mg2re in y we find 
y = 0.283 ± .069a; - 0.133 ± .007. The difference to the 
above correlation for central indices is barely significant. 

The quality of these correlations proves that both in- 
dices are essentially controlled by the same physical vari- 
ables, and leaves little room for the effects of diffuse dust 
upon the colours of E-galaxies. 

5. Conclusions 

New colour-radius relations have been derived for 36 E- 
type galaxies of the northern Local Supercluster, using 
UBVRI frames obtained with the 120cm telescope of the 
Observatoire de Haute-Provence. Four SAO, i.e. NGC3115, 
3607, 4551 and 5866 were also observed. 

We aimed to take advantage of the large field of the 
camera to observe the galaxies at larger radii than hitherto 
feasible, and thus improve the accuracy of colour gradi- 
ents. The availability of the series of aperture photometry 
in PP88 and PN94 for most of the sample was also consid- 
ered an asset towards a more coherent system of colours. 
It appears indeed that the colour calibrations are improved 
here compared to previous work, if this can be judged from 
the quality of correlations between zero point colours in 
various surveys (see Sect. 4.1) 

Two steps in the reduction procedure were thought 
significant in improving the quality of colour profiles: the 
first was the adjustment of the FWHM of the PSFs in a 
given colour set of 5 frames to the best of the five. This 
allowed us to get significant colours much closer to the 
galaxy center than otherwise feasible. The second was a 
careful " mapping" of the background of each frame, in or- 
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Table 9. Linear representation of colour against logr. Successive columns: NGC No.; Type; inner radius of calcu- 
lation; To outer radius of calculation; ro radius of colour evaluation (in log and "); U— B at ro and estimated standard 
error; Aj/s radial gradient; B— V at ro and estimated standard error; Asy radial gradient; V— R at ro and estimated 
standard error; Ayr radial gradient; V— I at ro and estimated standard error: Ay/ radial gradient; du Dust visibility 
index (Michard 1999) Notes: 2974: V-R and V-I not measurable; 3193: V-I not measurable. 
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der to lessen the background fluctuations remaining after 
the usual flat-fielding procedures. Both these precautions 
proved successful, and, as a result, the radial range of sat- 
isfactory colour measurements was greatly enlarged. Near 
the galaxy center, it proved feasible to obtain "central 
colours", i.e. colours integrated in the circle r = 3", in fair 
agreement with high resolution data (see Sect. 3.2.7.2 and 
Table 10). 



On the other hand, colours could be obtained at much 
lower surface brightness (or larger radii) than in previous 
work. Our colour data extend to 23.2 < /iy < 25.2, with a 
median value near jiy = 24.5 in all colours. According to 
the comparisons in Sect. 4.1, this is 1.5 to 2 magnitudes 
deeper than in previous work. " External colours" , refering 
to the level /xy = 24.5 whenever possible, are published 
for the first time (see Table 11), and may be useful to give 
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Table 10. Columns (1) to (5): integrated colours within 
a radius r = 3 ", i.e. UB^, BV3 , VR3 , VI3, UV3; (6): 
UVi,5 U— V colour at r = 1.5" calculated from Table 9 . 
(7) VR1.5 V-R colour at r = 1.5" 
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Table 11. "External colours" measured in the outer- 
most range of the data. Successive columns: fiy V surface 
brightness of measurement, a range of 0.5 magnitude cen- 
tered at the quoted beeing used. Columns (2), (3), (4) 
and (5): U— B, B— V, V— R, V— I respectively with esti- 
mated errors. These are the same in V— R as in B— V. 
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some indications about stellar populations at the outskirts 

of E-galaxies. 

On the other hand, the " red halo" effect of the camera 
was found to give enormous errors in V— I colours and gra- 
dients. These were corrected by a rigorous technique, and 
results in agreement with "classical" data were obtained. 
Considering the V— I gradients, one is not happy however 
to introduce in their evaluation, corrections larger than 
the quantity to be measured! Besides this specific prob- 



lem of the red halo of thinned CCD, the far wings of the 
PSFs have been proven in a recent paper (RMOl) to have 
non negligible effects in the gradients of other colours, and 
also to vary with the age of mirror coatings. It is not im- 
possible that the U— B or U— V gradients given here are 
overestimated by 15-20%, although they are in excellent 
statistical agreement with the well-known work of Peletier 
et al. (1990). 
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Various colour gradients against logr for a given ob- 
ject are well correlated, generally better than in previous 
work (see statistics in Sect. 4.1), which is interpreted as 
due to smaller measuring errors, notably in U— B. These 
improvements in accuracy did not bring out any obvious 
correlation between gradients and other galaxy properties. 
A few galaxies have exceptionally steep colour gradients 
(nearly at 2a) without sharing other properties. 

Colour-colour relations can be built from the present 
data for several locations in galaxies, such as near cen- 
ter, various fractions of the effective radius Vg, or the 
"outermost" measured range around = 24. All these 
diagrams overlap to form a single stripe with moderate 
scatter (except for one rather obvious calibration error?). 
These might prove useful to test theories of old stellar 
populations and of their host galaxies. Colour-colour di- 
agrams based upon integrated colours have already been 
used for this purpose (Worthey, 1994). 

The U— B or U— V colours correlate very well with the 
Mg2 index, both near the galaxy center and at the effec- 
tive radius rg- This seems to rule out any large influence 
of diffuse dust in the colours and colour gradients in E- 
galaxies. This was considered likely by Witt et al. (1992) 
and discussed by Wise and Silva (1996) with inconclusive 
results. Previous arguments against such an influence were 
presented in RMOO: they were based upon the relative av- 
erage values of the gradients in various colours, and are 
reinforced in the present work, since the mean gradients 
are nearly unchanged, and their errors lessened. 

It is well known that, for single-burst stellar popula- 
tions, colours and line indices depend both on the metallic- 
ity and on the age of the system (Worthey 1994; Borges et 
al. 1995). However, E-galaxies are constituted by a popu- 
lation mix, having age and metallicity distributions which 
reflect their star formation histories. Therefore a colour- 
metallicity calibration requires the use of models able to 
provide those distributions and, consequently, the inte- 
grated colours along the galaxy lifetime. Such a calibration 
will be presented in Paper II. 
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